Coulomb crystallization in two- component plasmas 
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The analysis of Coulomb crystallization is extended from one-component to two-component plas- 
mas. Critical parameters for the existence of Coulomb crystals are derived for both classical and 
quantum crystals. In the latter case, a critical mass ratio of the two charged components is found 
which is of the order of 80. Thus, holes in semiconductors with sufficiently flat valence bands are 
predicted to spontaneously order into a regular lattice. Such hole crystals are intimately related to 
ion Coulomb crystals in white dwarf and neutron stars as well as to ion crystals produced in the 
laboratory. A unified phase diagram of two-component Coulomb crystals is presented and is verified 
by first-principle computer simulations. 

PACS numbers: 52.27.-h, 52.25.-b 



Crystallization is one of the most fundamental many- 
particle phenomena in charged particle systems. After 
the prediction of a highly correlated state of the electron 
gas - the electron Wigner crystal [l| - there has been an 
active search for this phenomenon in many fields. Crys- 
tallization of electrons was observed on the surface of 
helium droplets [2| and is predicted to occur in semi- 
conductor quantum dots 3]. Moreover, crystals of ions 
have been observed in traps 4] and storage rings [B| , and 
are expected to occur in layered systems 0. The neces- 
sary condition for the existence of a crystal in these one- 
component plasmas (OCP) is that the mean Coulomb in- 
teraction energy, e 2 /f (f denotes the mean inter-particle 
distance), exceeds the mean kinetic energy (thermal en- 
ergy ^ksT or Fermi energy Ep in classical or quan- 
tum plasmas, respectively) by a factor T larger than T cr 
which, in a classical OCP is given by 175 @,L|. In a quan- 
tum OCP at zero temperature the coupling strength is 
measured by the Brueckner parameter, r s = f/a# (a# 
denotes the effective Bohr radius), with a critical value 
for crystallization of rf « 100 

The vast majority of Coulomb matter in the Universe, 
however, exists in the form of neutral plasmas, contain- 
ing (at least) two oppositely charged components (two- 
component plasma, TCP). Coulomb crystallization in a 
TCP has been observed as well, e.g. in colloidal and 
dusty plasmas ^| and it is predicted to be possible in 
laser-cooled expanding plasmas [ll| . The lattice of heavy 
particles is immersed into a structureless gas of the light 
component which does not affect the former. Besides 
these classical TCP crystals it is expected that in the in- 
terior of white dwarf stars and in the crust of neutron 
stars there exists an entirely different type of TCP crys- 
tals [l2|: Crystals of highly charged ions (e.g. fully ion- 
ized carbon, oxygen, iron) which are embedded into an 
extremely dense degenerate Fermi gas of electrons. No 
such quantum TCP crystals have been observed in the 
laboratory, despite early suggestions [l3|. It is an open 
question what classical and quantum TCP crystals (be- 



ing separated by 15 ... 20 orders of magnitude in density) 
have in common and if there exists an integrative phase 
diagram. 

This Letter aims to answer these questions. We show 
that, in fact, a common phase diagram of Coulomb crys- 
tals in a generic neutral TCP (consisting of electrons 
and point-like ions exists which is governed by five 
parameters - density and temperature (as in the OCP 
case) and, additionally, by the asymmetry of the heavy 
(h) and electron (e) components with respect to three 
fundamental properties: Mass, charge and temperature, 
M = m h /m e , Z = Z h /Z e and 6 = T e /T h . We show 
that classical TCP crystals require a critical charge ra- 
tio Z whereas quantum TCP crystals exist only if the 
mass ratio M exceeds a critical value of about 80. As a 
consequence, we predict the existence of quantum TCP 
crystals of protons in a dense hydrogen plasma and crys- 
tals of holes in semiconductors. 

Let us consider a locally neutral system of electrons 
and Z-fold charged heavy particles. The stationary state 
of the TCP is characterized by the dimensionless electron 
temperature T e — 3/c#X '/2Eb and mean inter-electron 
distance r se = f e /aB, where Eb denotes the e-h binding 
energy, and the dimensionless density is given by na B = 
3/(47ni). 

Classical crystal. For the existence of a Coulomb crys- 
tal in the presence of a classical gas of electrons, we first 
require that the heavy component is able to form a clas- 
sical OCP crystal, i.e. 



(i) 



and, secondly, that the electrons do not destroy that crys- 
tal, e.g., as a result of screening of the heavy particle 
interaction. However, the main obstacle for the crys- 
tal turns out to be the formation of e-h bound states 
(atoms, excitons etc.), because this drastically reduces 
the h-h correlation energy causing violation of condition 
(PP). Therefore, we require that no significant fraction of 
heavy particles is trapped in Coulomb bound states, for 
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which a conservative estimate is given by T e > 1, stating 
that the electrons have sufficiently high kinetic energy to 
escape the binding potential V e h- Making use of charge 
neutrality, n e = n^Z, and the definitions of T and T e we 
find from Eq. that the classical TCP Coulomb crystal 
exists between the temperatures and Tq 2 \ given by 



T^nl and T&\i 



36Z 2 / 3 1 



(2) 



and only if the ion charge exceeds a critical value, 
Z e cr (9,r se ) = (r cr r se /3B) 3/2 , which is independent of 
the mass ratio M. 

Quantum TCP crystal and critical mass ratio. In the 
presence of quantum electrons, the condition for crys- 
tallization of the heavy particles follows from the quan- 
tum OCP result, r s h > r c s r '. To be specific, we will con- 
centrate on hydrogen-like Coulomb bound states, where 

E b = 4^72^' and a B = £^ [er and m r are 
the background dielectric constant and the reduced mass 
m~ x = m^il + M)]. Since r c s r refers to the critical in- 
terparticle distance in units of the hydrogenic [8| , we 
first transform to the relevant effective Bohr radius as- 
Further, we eliminate r s h by expressing fu by f e using 
again charge neutrality. Then the OCP result can be 
rewritten as 



Z 4 / 3 (M + l)r se >r c s 



(3) 



This crystal of heavy particles will survive in the presence 
of electrons only if, as in the classical case, bound states 
are unstable which, at zero temperature, occurs due to 
pressure ionization at densities above the Mott density, 
i.e. r se < r^ ott ^1.2 (see below). With increasing tem- 
perature, ionization becomes possible at lower density 
which we indicate by a monotonically decreasing function 
l/r^ ott (T e ) which vanishes when T e — > 1 because there 
thermal ionization prevails. Combination of ((3J) with the 
existence of pressure ionization allows us to eliminate r se , 
yielding the criterion for the existence of a TCP crystal 
in the presence of a neutralizing background of quantum 
electrons as 



In absolute units, the density interval scales as Z 3 and the 
critical temperature as M • Z 4 • B. The critical mass ratio 
for singly charged ions in an isothermal TCP (B = 1) 
equals 83 and decreases with increasing temperature (due 
to the lower Mott density) and with increasing Z. 

Examples (0 = 1). (i) For a crystal of C 6+ [0 8+ ] ions 
expected to exist in the interior of white dwarf stars, 
the minimum density is given by n { e\0) = 2 • 10 26 cm 
[6.6 • 10 26 cm- 3 ], and T* = 10 9 K [4.2 • 10 9 K] 15]. (ii) 
Hydrogen and helium are predicted to form crystals as 
well: A crystal of protons [a— particles] is stable between 



ne '(Uj = U.9 • 10 24 cm- 3 [7.5 • 10 24 cm- 3 ] and n^ 2) (0) 

10 28 cm -3 [5 . l 30 cm -3] ? and below T * = g g . W 4 K 

[4.2 • 10 6 if] and should be achievable in laboratory ex- 
periments with laser or ion beam compression techniques, 
(hi) Crystallization of holes in semiconductors (Z = 1) is 
predicted for materials with a hole to electron mass ratio 
M > M cr « 83 [l6| . This value is feasible in intermedi- 
ate valence semiconductors, such as Tm[Se,Te] [TjJ. For 
example, for M = 100 (using e r = 20) the parameters 



are (0) 
T* w 9K. 
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FIG. 1: Snapshots of an e-h plasma at T e = 0.096, and 
r se — 0.63. Clouds of blue (yellow) dots mark the fully delo- 
calized electrons with spin up (down), clouds of red (pink) 
dots denote holes with spin up (down). Top left (right): 
M = 5 (50), bottom left (right): M = 100 (800). 



M > M cr (Z,T e ) 



Z 4/3 r Mott( Te ) 

which exists in a finite electron density range between 
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n (1) (T e ) 
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(M + 1)/(M C 



, n^(T e ) = n^(T e )K 3 , 
■1), (5) 



and below a critical temperature T*. The density limits 
follow from the Mott criterion and from Eq. (0 , whereas 
T* is estimated Q by the crossing point of the classical 
and quantum asymptotics of an OCP crystal (QJ and © 
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Z 2 0(M + 1) 



(6) 



Numerical verification. Of course, the boundary of the 
TCP crystal phase at low densities contains some uncer- 
tainty, due to the very complex transition from an atomic 
system to a Coulomb crystal of heavy particles embed- 
ded into delocalized electrons. This transition which ex- 
tends over a finite density interval may involve liquid- 
like behavior, clusters and, at low temperature, phase 
separation, an analysis of which is beyond the present 
study. We estimate that these effects give rise to an un- 
certainty of the minimum density (Mott density), 
of the order of 30%. Further, the error of rf is about 
20% Q , thus the critical parameters carry an uncertainty 
of about 50%. For particular systems, more accurate 
predictions are possible if the Mott parameter r^ ott is 
known, e.g. from computer simulations. Note that the 
complex processes of interest pose an extreme challenge 
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to the simulations: They must self-consistently include 
the full Coulomb interactions, e-h bound state formation 
in the presence of a surrounding plasma, pressure ion- 
ization and the quantum and spin properties of the light 
and heavy species. 

We therefore have performed extensive direct fermionic 
path integral Monte Carlo (PIMC) simulations of 
electron-hole plasmas which are based on our previous 
results for dense hydrogen-helium plasmas p^L e-h plas- 
mas JJj and electron Wigner crystallization While 
the so-called sign problem prohibits PIMC simulations 
of the ground state of a fermion system, here we restrict 
ourselves to temperatures at the upper boundary of the 
hole crystal phase, i.e. T e = 0.06 . . .0.2. Studying mass 
ratios in the range of M = 1 . . . 2000 and densities corre- 
sponding to r se = 0.6 ... 13 the simulations cover a large 
variety of Coulomb systems - from positronium, over typ- 
ical semiconductors to hydrogen. We start with the case 
of low densities (large r se ) to determine the Mott density 
n^. Here the TCP consists of excitons and biexcitons, 
and we found [2^j that for r se < 1.2 less than 10% of 
the electrons and holes are bound, approving the choice 



of 



^Mott 



made above. Thus crystallization should become 



possible. This is confirmed by our simulations, see Figs.Q] 
and showing results for r se = 0.63 and different mass 
ratios M. Fig. [I] displays snapshots of the e-h system in 
the simulation box. In all figure parts, the electrons form 
a nearly homogeneous Fermi gas - individual electrons 
penetrate each other, extending far beyond the main sim- 
ulation cell (shown by the grey grid lines). At the same 
time, the hole arrangement changes dramatically. While, 
at M = 5, the holes are in a gas-like state (similar to the 
electrons), at M = 50 their structure resembles a liquid 
and, at M = 100 and 800, they are periodically arranged 
in space. Thus, between M = 100 and 50 the holes crys- 
tallize. The figure also clearly shows the mechanism of 
this quantum melting process: With reduction of M the 
individual hole wave functions grow continuously until, 
at M = M cr , they exceed a critical size which allows for 
tunneling between lattice sites giving rise to hole der- 
ealization, i.e. crystal melting. Vice versa, increase of 
M reduces the spatial extension of each hole which sta- 
bilizes the crystal (at M = 800, they shrink to a dot, see 
Fig.DJ. 

The crystallization transition is further supported by 
the behavior of the pair distribution functions, g a 6, shown 
in Fig. [21 At M < 50 (upper two panels) the h-h func- 
tions (black curves) has only a single peak like in a liquid. 
However, for M = 100, ghh exhibits periodic oscillations 
with a deep first minimum, typical for a crystal. The 
crystal exists only at low temperature, an increase of T 
by a factor of 2 (from the lower right figure to the lower 
left) causes thermal melting. Further, the e-e and e-h pair 
distributions allow us to understand the behavior of the 
electrons: In the hole crystal phase the electrons exhibit 
periodic density modulations indicating the formation of 
Bloch waves (band structure) with increased concentra- 
tion at distances smaller than 0.5a#. Finally, we com- 



puted the relative distance fluctuations of the holes as a 
function of mass ratio at T e = 0.64 and r se = 0.63. They 
show a strong increase around M = 80 . . . 100 typical for 
a solid-liquid transition (not displayed). 
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FIG. 2: (Color online) e-e (blue, dash-dots), h-h (black, full 
line) and e-h (red, dotted line) pair distribution functions for 
different mass ratios and temperatures. 

Phase diagram. We now construct the generic phase 
diagram of the TCP which applies to all the different 
Coulomb crystals, see Fig. 03 Consider first the case 
of a hole crystal in semiconductors which is embedded 
into a dense Fermi gas of electrons. The holes behave 
classically above the black dotted line and quantum- 
mechanically below (this line is given by n^A| = 1, where 
Ah = h/ \/2irmhkBTh is the hole deBroglie wave length). 
The e-h bound state phase is shown in the left part and 
contains excitons and biexcitons and, eventually at low 
temperature, a Bose condensate, superfluid or an exci- 
tonic insulator [Tjj. On leaving this phase across its 
boundary (given by the blue line r se (T e ) = r^ ott (T e )) 
the fraction of bound states rapidly vanishes in favor 
of unbound e-h pairs with the holes showing liquid-like 
behavior. Upon further compression (at temperatures 
T e < T*) the hole liquid crystallizes, provided M > M cr . 
At the density quantum melting of the crystal is 
observed (vertical dashed green line). The entire hole 
Coulomb crystal phase for M=200 is marked by the full 
red line in Fig. [31 

Now, how is the hole crystal in semiconductors related 
to the Coulomb crystals in classical and astrophysical 
plasmas mentioned in the introduction? To answer this 
question we investigate the dependence of the stability 
of the crystal phase on the three asymmetry parame- 
ters. When M is reduced the crystal phase shrinks (see 
the red dashed line corresponding to M = 100) until 
for M = M cr it vanishes. Vice versa, when M becomes 
larger, both maximum density and temperature at which 
crystallization is possible, increase according to Eqs. (0 
E}. While in semiconductors in quasi-equilibrium M is 
the only parameter which can be varied from one ma- 
terial to another, the diversity of ionic plasmas, on the 
other hand, offers additionally control of the charge and 
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FIG. 3: (Color online) Phase diagram of a two- component 
plasma in the plane of dimensionless electron temperature 
T e and density parameter l/r se . The blue line indicates 
the boundary of the Coulomb bound state phase given by 
r^ ott (T e ). Above (below) the dotted black line holes are clas- 
sical (degenerate). The red full (dashed) line is the boundary 
of the hole crystal for 6 = Z = 1 and M = 200 (M 100) 
with the asymptotics given by Eqs. (|1I5|) (green dashed lines). 

temperature ratios Z and O in very broad ranges. By 
increasing M, Z and 6 the crystal phase extends further 
towards higher density and temperature covering an ever 
increasing part of the temperature-density plane (with 
the exception of the bound state phase). Eventually this 
crystal phase will overlap with the known classical and 
astrophysical Coulomb crystals at low and high densi- 
ties, respectively. Thus, indeed, the phase diagram in 
Fig. 01 applies to all Coulomb crystals in two-component 
plasmas of electrons and point-like ions, independently 
of their physical origin. Of course, in different systems 
specific additional factors may exist. For example, phase 
separation or non-Coulombic bound states will modify 
the boundary of the bound state phase and the value 



of whereas band structure effects or extended ionic 
cores can modify the high-density behavior 15]. Finally, 
while our analytical results for the crystal phase are ob- 
tained neglecting e-h correlations (which is justified by 
the large mass ratio M) , the simulations indicate that the 
electrons have a stabilizing influence, increasing the max- 
imum temperature beyond T*. These effects are beyond 
the present investigation and will be discussed elsewhere. 

In summary, one of the most fundamental collective 
properties of matter - Coulomb crystallization has been 
extended from the OCP model to the case of neutral 
systems with two oppositely charged components. Our 
analysis provides a common view and general quantita- 
tive bounds on the critical parameters for the existence 
of Coulomb crystals in a large variety of TCP, including 
dwarf stars, laser-cooled expanding plasmas, dusty plas- 
mas and semiconductors. The critical parameters depend 
on the combined mass, charge and temperature asym- 
metry between the heavy and light component. Crys- 
tallization of protons and of holes in semiconductors is 
predicted. Hole crystals should exist in materials with 
a mass ratio of about 80 21] and might be observable 
in rare earth semiconductors in neutron scattering ex- 
periments. These hole crystals could serve as a valuable 
testing ground for quantum TCP crystals in general and 
for ion crystals in exotic stellar objects in particular. 
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